Translational modeling-based evidence for enhanced efficacy of standard-of-care drugs in combination with anti-microRNA-155 in non-small-cell lung cancer.

Elevated microRNA-155 (miR-155) expression in non-small-cell lung cancer (NSCLC) promotes cisplatin resistance and negatively impacts treatment outcomes. However, miR-155 can also boost anti-tumor immunity by suppressing PD-L1 expression. We developed a multiscale mechanistic model, calibrated with in vivo data and then extrapolated to humans, to investigate the therapeutic effects of nanoparticle-delivered anti-miR-155 in NSCLC, alone or in combination with standard-of-care drugs. Model simulations and analyses of the clinical scenario revealed that monotherapy with anti-miR-155 at a dose of 2.5 mg/kg administered once every three weeks has substantial anti-cancer activity. It led to a median progression-free survival (PFS) of 6.7 months, which compared favorably to cisplatin and immune checkpoint inhibitors. Further, we explored the combinations of anti-miR-155 with standard-of-care drugs, and found strongly synergistic two- and three-drug combinations. A three-drug combination of anti-miR-155, cisplatin, and pembrolizumab resulted in a median PFS of 13.1 months, while a two-drug combination of anti-miR-155 and cisplatin resulted in a median PFS of 11.3 months, which emerged as a more practical option due to its simple design and cost-effectiveness. Our analyses also provided valuable insights into unfavorable dose ratios for drug combinations, highlighting the need for optimizing dose regimen to prevent antagonistic effects. Thus, this work bridges the gap between preclinical development and clinical translation of anti-miR-155 and unravels the potential of anti-miR-155 combination therapies in NSCLC.

Hybridizing mechanistic mathematical modeling with deep learning methods to predict individual cancer patient survival after immune checkpoint inhibitor therapy.

We present a study where predictive mechanistic modeling is used in combination with deep learning methods to predict individual patient survival probabilities under immune checkpoint inhibitor (ICI) therapy. This hybrid approach enables prediction based on both measures that are calculable from mechanistic models (but may not be directly measurable in the clinic) and easily measurable quantities or characteristics (that are not always readily incorporated into predictive mechanistic models). The mechanistic model we have applied here can predict tumor response from CT or MRI imaging based on key mechanisms underlying checkpoint inhibitor therapy, and in the present work, its parameters were combined with readily-available clinical measures from 93 patients into a hybrid training set for a deep learning time-to-event predictive model. Analysis revealed that training an artificial neural network with both mechanistic modeling-derived and clinical measures achieved higher per-patient predictive accuracy based on event-time concordance, Brier score, and negative binomial log-likelihood-based criteria than when only mechanistic model-derived values or only clinical data were used. Feature importance analysis revealed that both clinical and model-derived parameters play prominent roles in neural network decision making, and in increasing prediction accuracy, further supporting the advantage of our hybrid approach. We anticipate that many existing mechanistic models may be hybridized with deep learning methods in a similar manner to improve predictive accuracy through addition of additional data that may not be readily implemented in mechanistic descriptions.

Mechanistic modeling for optimal design of dissolvable microneedle-based patches for transdermal drug delivery.

Polymeric microneedle (MN)-based patches are an efficient, non-invasive, and painless means of drug delivery through the skin to systemic circulation. The design of these MN-based patches can be customized for various drug delivery applications, particularly modified release of drugs. In this study, we developed a mathematical model of drug delivery via MN-based patches to study the effect of patch design properties on drug delivery kinetics and systemic pharmacokinetics (PK). We calibrated the model against two representative formulations: a rapid release patch of naloxone and a sustained-release patch of levonorgestrel. The model was then applied to assess the relative significance of model parameters in governing systemic PK of drugs and obtain optimal design parameters to achieve therapeutically meaningful drug levels in a clinical setting. We identified the importance of drug loading fraction, MN base radius, and MN height as the key control parameters responsible for drug PK.Clinical Relevance- Through the application of modeling and simulation, we can improve drug delivery from MN-based patches by identifying optimal design parameters to support the clinical translation of these novel drug delivery systems.

Hybrid Cellular Automata Modeling Reveals the Effects of Glucose Gradients on Tumour Spheroid Growth.

PURPOSE: In recent years, mathematical models have become instrumental in cancer research, offering insights into tumor growth dynamics, and guiding the development of pharmacological strategies. These models, encompassing diverse biological and physical processes, are increasingly used in clinical settings, showing remarkable predictive precision for individual patient outcomes and therapeutic responses. METHODS: Motivated by these advancements, our study introduces an innovative in silico model for simulating tumor growth and invasiveness. The automated hybrid cell emulates critical tumor cell characteristics, including rapid proliferation, heightened motility, reduced cell adhesion, and increased responsiveness to chemotactic signals. This model explores the potential evolution of 3D tumor spheroids by manipulating biological parameters and microenvironment factors, focusing on nutrient availability. RESULTS: Our comprehensive global and local sensitivity analysis reveals that tumor growth primarily depends on cell duplication speed and cell-to-cell adhesion, rather than external chemical gradients. Conversely, tumor invasiveness is predominantly driven by chemotaxis. These insights illuminate tumor development mechanisms, providing vital guidance for effective strategies against tumor progression. Our proposed model is a valuable tool for advancing cancer biology research and exploring potential therapeutic interventions.

Mechanistic modeling of anti-microRNA-155 therapy combinations in lung cancer.

To improve treatment outcomes in non-small cell lung cancer (NSCLC), it is crucial to identify treatment strategies with the potential to exhibit drug synergism. This can lower the required effective dose, reducing exposure to drugs and associated toxicities, while improving treatment efficacy. In previous studies, drugs targeting the microRNA-155 or PD-L1 have been promising in restraining NSCLC tumor growth. We have developed a mathematical model that simulates the in vivo pharmacokinetics and pharmacodynamics of the novel nanoparticle-delivered anti-microRNA-155 for potential use with standard-of-care drug atezolizumab for NSCLC. Through modeling and simulation, we identified possible drug synergism between the two drugs that holds promise to improve tumor response at reduced drug exposure.Clinical Relevance-Identifying the possibility of drug synergism for an anti-microRNA-155 based nanotherapeutic with standard-of-care immunotherapy to improve lung cancer treatment outcomes.

Exploring Cell Migration Mechanisms in Cancer: From Wound Healing Assays to Cellular Automata Models.

PURPOSE: Cell migration is a critical driver of metastatic tumor spread, contributing significantly to cancer-related mortality. Yet, our understanding of the underlying mechanisms remains incomplete. METHODS: In this study, a wound healing assay was employed to investigate cancer cell migratory behavior, with the aim of utilizing migration as a biomarker for invasiveness. To gain a comprehensive understanding of this complex system, we developed a computational model based on cellular automata (CA) and rigorously calibrated and validated it using in vitro data, including both tumoral and non-tumoral cell lines. Harnessing this CA-based framework, extensive numerical experiments were conducted and supported by local and global sensitivity analyses in order to identify the key biological parameters governing this process. RESULTS: Our analyses led to the formulation of a power law equation derived from just a few input parameters that accurately describes the governing mechanism of wound healing. This groundbreaking research provides a powerful tool for the pharmaceutical industry. In fact, this approach proves invaluable for the discovery of novel compounds aimed at disrupting cell migration, assessing the efficacy of prospective drugs designed to impede cancer invasion, and evaluating the immune system's responses.

Ceramide as an endothelial cell surface receptor and a lung-specific lipid vascular target for circulating ligands.

The vascular endothelium from individual organs is functionally specialized, and it displays a unique set of accessible molecular targets. These serve as endothelial cell receptors to affinity ligands. To date, all identified vascular receptors have been proteins. Here, we show that an endothelial lung-homing peptide (CGSPGWVRC) interacts with C16-ceramide, a bioactive sphingolipid that mediates several biological functions. Upon binding to cell surfaces, CGSPGWVRC triggers ceramide-rich platform formation, activates acid sphingomyelinase and ceramide production, without the associated downstream apoptotic signaling. We also show that the lung selectivity of CGSPGWVRC homing peptide is dependent on ceramide production in vivo. Finally, we demonstrate two potential applications for this lipid vascular targeting system: i) as a bioinorganic hydrogel for pulmonary imaging and ii) as a ligand-directed lung immunization tool against COVID-19. Thus, C16-ceramide is a unique example of a lipid-based receptor system in the lung vascular endothelium targeted in vivo by circulating ligands such as CGSPGWVRC.

Prospective Subunit Nanovaccine against Mycobacterium tuberculosis Infection─Cubosome Lipid Nanocarriers of Cord Factor, Trehalose 6,6' Dimycolate.

An improved vaccine is urgently needed to replace the now more than 100-year-old Bacillus Calmette-Guérin (BCG) vaccine against tuberculosis (TB) disease, which represents a significant burden on global public health. Mycolic acid, or cord factor trehalose 6,6' dimycolate (TDM), a lipid component abundant in the cell wall of the pathogen Mycobacterium tuberculosis (MTB), has been shown to have strong immunostimulatory activity but remains underexplored due to its high toxicity and poor solubility. Herein, we employed a novel strategy to encapsulate TDM within a cubosome lipid nanocarrier as a potential subunit nanovaccine candidate against TB. This strategy not only increased the solubility and reduced the toxicity of TDM but also elicited a protective immune response to control MTB growth in macrophages. Both pre-treatment and concurrent treatment of the TDM encapsulated in lipid monoolein (MO) cubosomes (MO-TDM) (1 mol %) induced a strong proinflammatory cytokine response in MTB-infected macrophages, due to epigenetic changes at the promoters of tumor necrosis factor alpha (TNF-α) and interleukin 6 (IL-6) in comparison to the untreated control. Furthermore, treatment with MO-TDM (1 mol %) cubosomes significantly improved antigen processing and presentation capabilities of MTB-infected macrophages to CD4 T cells. The ability of MO-TDM (1 mol %) cubosomes to induce a robust innate and adaptive response in vitro was further supported by a mathematical modeling study predicting the vaccine efficacy in vivo. Overall, these results indicate a strong immunostimulatory effect of TDM when delivered through the lipid nanocarrier, suggesting its potential as a novel TB vaccine.

A modeling-based approach to optimize COVID-19 vaccine dosing schedules for improved protection.

While the development of different vaccines slowed the dissemination of SARS-CoV-2, the occurrence of breakthrough infections has continued to fuel the COVID-19 pandemic. To secure at least partial protection in the majority of the population through 1 dose of a COVID-19 vaccine, delayed administration of boosters has been implemented in many countries. However, waning immunity and emergence of new variants of SARS-CoV-2 suggest that such measures may induce breakthrough infections due to intermittent lapses in protection. Optimizing vaccine dosing schedules to ensure prolonged continuity in protection could thus help control the pandemic. We developed a mechanistic model of immune response to vaccines as an in silico tool for dosing schedule optimization. The model was calibrated with clinical data sets of acquired immunity to COVID-19 mRNA vaccines in healthy and immunocompromised participants and showed robust validation by accurately predicting neutralizing antibody kinetics in response to multiple doses of COVID-19 mRNA vaccines. Importantly, by estimating population vulnerability to breakthrough infections, we predicted tailored vaccination dosing schedules to minimize breakthrough infections, especially for immunocompromised individuals. We identified that the optimal vaccination schedules vary from CDC-recommended dosing, suggesting that the model is a valuable tool to optimize vaccine efficacy outcomes during future outbreaks.

Mechanistic modeling-guided optimization of microneedle-based skin patch for rapid transdermal delivery of naloxone for opioid overdose treatment.

Naloxone, an FDA-approved opioid inhibitor, used to reverse opioid overdose complications has up till date faced challenges associated with its delivery. Limitations include the use of invasive delivery forms and the need for frequent redosing due to its short half-life. The goal of the current study was to design a transdermal rapidly dissolving polymeric microneedle (MN) patch with delivery and pharmacokinetic properties comparable to that seen with the commercially available NAL products, eliminating their delivery limitations. Patches of varying dimensions (500 µm; 100 array,800 µm; 100array, and 600 µm; 225 array) were fabricated to evaluate the effect of increasing MN length, and density (no. of needles/unit area) on drug release. Drug dose in each of these patches was 17.89 ± 0.23 mg, 17.2 ± 0.77 mg, and 17.8 ± 1.01 mg, respectively. Furthermore, the insertion efficiency of each of the MN patches was 94 ± 4.8%, 90.6 ± 1.69%, and 96 ± 1.29%, respectively. Compared to passive permeation, a reduced lag time of about 5-15 min was observed with a significant drug flux of 15.09 ± 7.68 g[Formula: see text]/cm2/h seen in the first 1 h (p < 0.05) with the array of 100 needles (500 µm long). Over 24 h, a four and ten-fold increase in permeation was seen with the longer length and larger density MN patch, respectively, when compared to the 500 µm (100 array) patch. Model simulations and analyses revealed the significance of needle base diameter and needle count in improving systemic pharmacokinetics of NAL.

Immune-checkpoint inhibitor therapy response evaluation using oncophysics-based mathematical models.

The field of oncology has transformed with the advent of immunotherapies. The standard of care for multiple cancers now includes novel drugs that target key checkpoints that function to modulate immune responses, enabling the patient's immune system to elicit an effective anti-tumor response. While these immune-based approaches can have dramatic effects in terms of significantly reducing tumor burden and prolonging survival for patients, the therapeutic approach remains active only in a minority of patients and is often not durable. Multiple biological investigations have identified key markers that predict response to the most common form of immunotherapy-immune checkpoint inhibitors (ICI). These biomarkers help enrich patients for ICI but are not 100% predictive. Understanding the complex interactions of these biomarkers with other pathways and factors that lead to ICI resistance remains a major goal. Principles of oncophysics-the idea that cancer can be described as a multiscale physical aberration-have shown promise in recent years in terms of capturing the essence of the complexities of ICI interactions. Here, we review the biological knowledge of mechanisms of ICI action and how these are incorporated into modern oncophysics-based mathematical models. Building on the success of oncophysics-based mathematical models may help to discover new, rational methods to engineer immunotherapy for patients in the future. This article is categorized under: Therapeutic Approaches and Drug Discovery > Nanomedicine for Oncologic Disease.

A modeling-based approach to optimize COVID-19 vaccine dosing schedules for improved protection.

While the development of different vaccines has slowed the dissemination of SARS-CoV-2, the occurrence of breakthrough infections continues to fuel the pandemic. As a strategy to secure at least partial protection, with a single dose of a given COVID-19 vaccine to maximum possible fraction of the population, delayed administration of subsequent doses (or boosters) has been implemented in many countries. However, waning immunity and emergence of new variants of SARS-CoV-2 suggest that such measures may jeopardize the attainment of herd immunity due to intermittent lapses in protection. Optimizing vaccine dosing schedules could thus make the difference between periodic occurrence of breakthrough infections or effective control of the pandemic. To this end, we have developed a mechanistic mathematical model of adaptive immune response to vaccines and demonstrated its applicability to COVID-19 mRNA vaccines as a proof-of-concept for future outbreaks. The model was thoroughly calibrated against multiple clinical datasets involving immune response to SARS-CoV-2 infection and mRNA vaccines in healthy and immunocompromised subjects (cancer patients undergoing therapy); the model showed robust clinical validation by accurately predicting neutralizing antibody kinetics, a correlate of vaccine-induced protection, in response to multiple doses of mRNA vaccines. Importantly, we estimated population vulnerability to breakthrough infections and predicted tailored vaccination dosing schedules to maximize protection and thus minimize breakthrough infections, based on the immune status of a sub-population. We have identified a critical waiting window for cancer patients (or, immunocompromised subjects) to allow recovery of the immune system (particularly CD4+ T-cells) for effective differentiation of B-cells to produce neutralizing antibodies and thus achieve optimal vaccine efficacy against variants of concern, especially between the first and second doses. Also, we have obtained optimized dosing schedules for subsequent doses in healthy and immunocompromised subjects, which vary from the CDC-recommended schedules, to minimize breakthrough infections. The developed modeling tool is based on generalized adaptive immune response to antigens and can thus be leveraged to guide vaccine dosing schedules during future outbreaks.

Zonal Intratumoral Delivery of Nanoparticles Guided by Surface Functionalization.

Delivery of small molecules and anticancer agents to malignant cells or specific regions within a tumor is limited by penetration depth and poor spatial drug distribution, hindering anticancer efficacy. Herein, we demonstrate control over gold nanoparticle (GNP) penetration and spatial distribution across solid tumors by administering GNPs with different surface chemistries at a constant injection rate via syringe pump. A key finding in this study is the discovery of different zone-specific accumulation patterns of intratumorally injected nanoparticles dependent on surface functionalization. Computed tomography (CT) imaging performed in vivo of C57BL/6 mice harboring Lewis lung carcinoma (LLC) tumors on their flank and gross visualization of excised tumors consistently revealed that intratumorally administered citrate-GNPs accumulate in particle clusters in central areas of the tumor, while GNPs functionalized with thiolated phosphothioethanol (PTE-GNPs) and thiolated polyethylene glycol (PEG-GNPs) regularly accumulate in the tumor periphery. Further, PEG functionalization resulted in larger tumoral surface coverage than PTE, reaching beyond the outer zone of the tumor mass and into the surrounding stroma. To understand the dissimilarities in spatiotemporal evolution across the different GNP surface chemistries, we modeled their intratumoral transport with reaction-diffusion equations. Our results suggest that GNP surface passivation affects nanoparticle reactivity with the tumor microenvironment, leading to differential transport behavior across tumor zones. The present study provides a mechanistic understanding of the factors affecting spatiotemporal distribution of nanoparticles in the tumor. Our proof of concept of zonal delivery within the tumor may prove useful for directing anticancer therapies to regions of biomarker overexpression.

Localization of drug biodistribution in a 3D-bioengineered subcutaneous neovascularized microenvironment.

Local immunomodulation has shown the potential to control the immune response in a site-specific manner for wound healing, cancer, allergy, and cell transplantation, thus abrogating adverse effects associated with systemic administration of immunotherapeutics. Localized immunomodulation requires confining the biodistribution of immunotherapeutics on-site for maximal immune control and minimal systemic drug exposure. To this end, we developed a 3D-printed subcutaneous implant termed 'NICHE', consisting of a bioengineered vascularized microenvironment enabled by sustained drug delivery on-site. The NICHE was designed as a platform technology for investigating local immunomodulation in the context of cell therapeutics and cancer vaccines. Here we studied the ability of the NICHE to localize the PK and biodistribution of different model immunomodulatory agents in vivo. For this, we first performed a mechanistic evaluation of the microenvironment generated within and surrounding the NICHE, with emphasis on the parameters related to molecular transport. Second, we longitudinally studied the biodistribution of ovalbumin, cytotoxic T lymphocyte-associated antigen-4-Ig (CTLA4Ig), and IgG delivered locally via NICHE over 30 days. Third, we used our findings to develop a physiologically-based pharmacokinetic (PBPK) model. Despite dense and mature vascularization within and surrounding the NICHE, we showed sustained orders of magnitude higher molecular drug concentrations within its microenvironment as compared to systemic circulation and major organs. Further, the PBPK model was able to recapitulate the biodistribution of the 3 molecules with high accuracy (r â€‹> â€‹0.98). Overall, the NICHE and the PBPK model represent an adaptable platform for the investigation of local immunomodulation strategies for a wide range of biomedical applications.

Genetic and Structural Analysis of SARS-CoV-2 Spike Protein for Universal Epitope Selection.

Evaluation of immunogenic epitopes for universal vaccine development in the face of ongoing SARS-CoV-2 evolution remains a challenge. Herein, we investigate the genetic and structural conservation of an immunogenically relevant epitope (C662-C671) of spike (S) protein across SARS-CoV-2 variants to determine its potential utility as a broad-spectrum vaccine candidate against coronavirus diseases. Comparative sequence analysis, structural assessment, and molecular dynamics simulations of C662-C671 epitope were performed. Mathematical tools were employed to determine its mutational cost. We found that the amino acid sequence of C662-C671 epitope is entirely conserved across the observed major variants of SARS-CoV-2 in addition to SARS-CoV. Its conformation and accessibility are predicted to be conserved, even in the highly mutated Omicron variant. Costly mutational rate in the context of energy expenditure in genome replication and translation can explain this strict conservation. These observations may herald an approach to developing vaccine candidates for universal protection against emergent variants of coronavirus.

Dedifferentiation-mediated stem cell niche maintenance in early-stage ductal carcinoma in situ progression: insights from a multiscale modeling study.

We present a multiscale agent-based model of ductal carcinoma in situ (DCIS) to study how key phenotypic and signaling pathways are involved in the early stages of disease progression. The model includes a phenotypic hierarchy, and key endocrine and paracrine signaling pathways, and simulates cancer ductal growth in a 3D lattice-free domain. In particular, by considering stochastic cell dedifferentiation plasticity, the model allows for study of how dedifferentiation to a more stem-like phenotype plays key roles in the maintenance of cancer stem cell populations and disease progression. Through extensive parameter perturbation studies, we have quantified and ranked how DCIS is sensitive to perturbations in several key mechanisms that are instrumental to early disease development. Our studies reveal that long-term maintenance of multipotent stem-like cell niches within the tumor are dependent on cell dedifferentiation plasticity, and that disease progression will become arrested due to dilution of the multipotent stem-like population in the absence of dedifferentiation. We have identified dedifferentiation rates necessary to maintain biologically relevant multipotent cell populations, and also explored quantitative relationships between dedifferentiation rates and disease progression rates, which may potentially help to optimize the efficacy of emerging anti-cancer stem cell therapeutics.

Novel intranasal treatment for anxiety disorders using amiloride, an acid-sensing ion channel antagonist: Pharmacokinetic modeling and simulation.

OBJECTIVE: To develop a physiologically based pharmacokinetic (PBPK) model for amiloride, an acid-sensing ion channel (ASIC) antagonist, and to simulate its pharmacokinetics in plasma and the central nervous system following intranasal administration in a virtual human population. MATERIALS AND METHODS: We first developed a PBPK model of amiloride after oral administration and optimized the model using data from five clinical studies. Next, we added a nasal compartment to the amiloride oral PBPK model and parameterized using data from previous clinical studies. We simulated amiloride's pharmacokinetics in plasma, brain, and cerebrospinal fluid (CSF) after intranasal administration of amiloride at various doses in a virtual human population. RESULTS: The target amiloride concentration in the central nervous system required for maximal ASIC inhibition was achieved with a 75-mg intranasal amiloride dose. However, this finding is based on simulations performed using a mathematical model and needs to be further validated with appropriate clinical data. CONCLUSION: The nasal PBPK model of amiloride could be used to design future clinical studies and allow for successful clinical translation of intranasal amiloride formulation.

Translational Modeling Identifies Synergy between Nanoparticle-Delivered miRNA-22 and Standard-of-Care Drugs in Triple-Negative Breast Cancer.

PURPOSE: Downregulation of miRNA-22 in triple-negative breast cancer (TNBC) is associated with upregulation of eukaryotic elongation 2 factor kinase (eEF2K) protein, which regulates tumor growth, chemoresistance, and tumor immunosurveillance. Moreover, exogenous administration of miRNA-22, loaded in nanoparticles to prevent degradation and improve tumor delivery (termed miRNA-22 nanotherapy), to suppress eEF2K production has shown potential as an investigational therapeutic agent in vivo. METHODS: To evaluate the translational potential of miRNA-22 nanotherapy, we developed a multiscale mechanistic model, calibrated to published in vivo data and extrapolated to the human scale, to describe and quantify the pharmacokinetics and pharmacodynamics of miRNA-22 in virtual patient populations. RESULTS: Our analysis revealed the dose-response relationship, suggested optimal treatment frequency for miRNA-22 nanotherapy, and highlighted key determinants of therapy response, from which combination with immune checkpoint inhibitors was identified as a candidate strategy for improving treatment outcomes. More importantly, drug synergy was identified between miRNA-22 and standard-of-care drugs against TNBC, providing a basis for rational therapeutic combinations for improved response CONCLUSIONS: The present study highlights the translational potential of miRNA-22 nanotherapy for TNBC in combination with standard-of-care drugs.

Mathematical modeling of cancer immunotherapy for personalized clinical translation.

Encouraging advances are being made in cancer immunotherapy modeling, especially in the key areas of developing personalized treatment strategies based on individual patient parameters, predicting treatment outcomes and optimizing immunotherapy synergy when used in combination with other treatment approaches. Here we present a focused review of the most recent mathematical modeling work on cancer immunotherapy with a focus on clinical translatability. It can be seen that this field is transitioning from pure basic science to applications that can make impactful differences in patients' lives. We discuss how researchers are integrating experimental and clinical data to fully inform models so that they can be applied for clinical predictions, and present the challenges that remain to be overcome if widespread clinical adaptation is to be realized. Lastly, we discuss the most promising future applications and areas that are expected to be the focus of extensive upcoming modeling studies.

A Multiscale Model to Identify Limiting Factors in Nanoparticle-Based miRNA Delivery for Tumor Inhibition.

MicroRNA-based gene therapy for cancer treatment via nanoparticles (NPs) requires navigation of multiple physical and physiological barriers in order to efficiently deliver the miRNAs to the cancer cell cytoplasm. We here present a mathematical model to investigate the variability associated with tumor, NP, and miRNA characteristics, and identify the limiting factors in miRNA delivery to tumors. Through global parameter analysis, the miRNA release rate from NPs and NP degradability were found to have the most significant impact on cytosolic accumulation of miRNAs. These NP properties can be fine-tuned in order to optimize the delivery system for enhancing the efficacy of miRNA-based therapy.Clinical Relevance-Understanding the effect of nanoparticle, tumor, and miRNA characteristics in governing the efficacy of miRNA-based cancer therapy will support its clinical translation.